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ABSTRACT 

Lattice QCD simulations with staggered fermions rely on the "fourth- 
root trick." The validity of this trick has been proved for free staggered 
fermions using renormalization-group block transformations. I review the 
elements of the construction and discuss how it might be generalized to 
the interacting case. 
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1. The fourth-root trick. A staggered-fermion field x( x ) consists of a single fermionic 
degree of freedom per color per lattice site [I]. The sixteen variables (per color) 
residing in each 2 4 hypercube of the euclidean lattice can be re-grouped into four 
quark fields in the continuum limit. Thus, a single staggered field could be used to 
simulate the up, down, strange and charm quarks, all at once, provided an appropriate 
mass matrix is chosen [2] . For various technical reasons, however, this is not the way 
staggered-fermion simulations are done. 

The "rooted" Boltzmann weight used to generate the dynamical staggered-fermion 
configurations contains a factor of det 1 ^ 4 (D s + m) for each of the three light quarks, 
where D s is the massless, interacting, anti-hermitian, one- component, staggered- 
fermion operator. 1 The raison d'etre behind the fourth- root trick is simple [3J. In 
the continuum limit, this staggered-fermion determinant will describe four equal-mass 
quarks that interact only through the exchange of gluons; hence, the determinant's 
fourth root is expected to account for a single quark with the same mass. The diffi- 
cultly is that, for any finite lattice spacing a, the four quark species (or four "tastes," 
reserving the term "flavors" to the different-mass species) are entangled at short dis- 
tances. There is no simple way to represent the rooted Boltzmann weight as a path 
integral with a local action. This raises the question whether the rooted staggered 
theory is consistent with the rules of local quantum field theory, and whether its con- 
tinuum limit is in the same universality class as QCD. The unprecedented accuracy 
of staggered-fermions simulations has made this issue all the more imminent [I]. 

In the free theory, the re-grouping alluded to in the first paragraph is explicitly 
given by [5] 

M*)= E (7[ 1 7 2 r2 73 3 7?)a 4 x(25 + r). (1) 

J>=0,1 

Here a and i are respectively Dirac and taste indices, that both run from 1 to 4. The 
taste-representation variables ip a i live on a coarse lattice, with spacing 2a, whose sites 
are labeled by the coordinates x. Each site of the original, fine lattice has a unique 
representation as x = 2x + r, where r M = 0, 1. In the taste basis the free (massless) 
staggered action is So = il>D ip, where 

D = a' 1 ([7m ® sin(p M a) + [75 <g> r 5 r /Lt ](l - cos(p M a)) , (2) 
n 

in momentum space. The two sets of Dirac matrices, 7^ and r M = 7* = 7J, act on 
the Dirac and taste indices respectively. I is the identity matrix in taste space. The 
kinetic term — the first term on the right-hand side — is the naive discretization of the 
continuum p. The second term lifts the doublers at the p 7^ corners of the Brillouin 
zone (of the coarse lattice). We will refer to it as a skewed — because of the extra r M 
matrices — Wilson term. 

As a start, there are two complementary ways to convince oneself that the fourth 
root "has to" work in the free theory. First, one can look for the four-fold degener- 
acy expected from the taste structure. In momentum space, the continuum limit of 
a free theory corresponds to \p\a — ► 0. The kinetic term starts off as f 1 + 0(a 2 p 3 ), 

1 In practice the one usually sets m u = nid = mii g ht and simulates det 1 ^ 2 (D s +mii g ht) det 1 ^ 4 (Z) s + 

^strange) ■ 
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whereas the skewed Wilson term starts off as 0(ap 2 ); thus the skewed Wilson term 
drops out, and the expected four-fold degeneracy is obtained. Equivalently, one can 
examine the behavior of the propagator Go = Dq 1 at large space-time separations. 
The momentum-space representation of Go is easily worked out. The only singularity, 
occurring at p — > 0, has the same structure as in the continuum. The taste-breaking 
terms, originating in the skewed Wilson term, do not produce a singularity in the 
propagator. As a result P, the taste-breaking part of the coordinate-space free prop- 
agator vanishes exponentially with the separation, with an 0(1) decay rate in lattice 
units. We must keep in mind, though, that for momenta or separations which are 
0(1) in lattice units, the skewed Wilson term does spoil the diagonal taste structure. 
Therefore we should not expect that merely taking the fourth root of the product of 
all the eigenvalues, small and large alike, will yield a sensible lattice operator. 

Next let us examine one-loop perturbation theory, taking the vacuum polarization 
(p) as an example. We recall that the interacting theory is defined using the one- 
component formulation, and that all of its symmetries are needed to ensure that mass 
terms will always renormalize multiplicatively |2 . How to construct an interacting 
theory in the taste basis is an issue that we discuss later. For massless staggered 
fermions one obtains 

n Mi /(p) = (Sf, u p 2 - p M p I/ )(4clog(po) + contact terms) . (3) 

The logarithmic term is universal. The coefficient c is the same as one would obtain for 
a single quark in a continuum calculation, or for that matter, in a lattice calculation 
using any one-flavor Dirac operator, such as a Wilson fermion. The non-analytic, 
long-distance part of the vacuum polarization of staggered fermions will therefore be 
reproduced by (say) four Wilson fermions. Taking the fourth root would replace these 
four Wilson fermions by one, implying that the fourth root of the long-distance part 
of the determinant is a sensible quantity. 

Like in any perturbative lattice calculation, one-loop diagrams with staggered 
fermions inside the loop also contain contact terms, which are all local and gauge 
invariant. The first of those will correspond to (a lattice discretization of) F 2 , and 
has the effect of a finite renormalization of the coupling constant. Other terms corre- 
spond to irrelevant operators. Similar statements apply to a calculation with Wilson 
fermions. However, the contact terms are not universal, and the staggered-fermions 
contact terms are not equal to four times the Wilson-fermions contact terms. But 
then, this short- distance discrepancy is not a problem: In lattice perturbation theory, 
the fourth root of the staggered determinant will be almost equal to a Wilson-fermion 
determinant; and, if we wish to, we can always make up for the difference in the con- 
tact terms by introducing innocuous, local modifications to the gauge-field action. 

The fourth-root trick is also valid within the framework of low-energy effective 
lagrangians that capture the infra-red limit of QCD [7] . The remaining challenge is to 
promote the arguments to a fully non-perturbative setup. The lesson is that we will 
need a device that can separate out long-distance from short-distance physics, thus 
allowing us to treat differently the two parts of the determinant. This is precisely 
what Renormalization-Group (RG) blocking was designed for jHJ E] . 
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2. Renormalization- group transformations. Our first task is to establish the validity 
of of the fourth-root trick in the free theory, using RG block transformations PU| ■ A 
single RG blocking step works as follows jHj: 

Z = J dipdifj exp(—if)D ip) 

= J dipdipdqdq exp (^—ipD ip — a(q — if>Q^){q — Qif>)^ 

= det(G^) J dqdq exp(-g£>ig) . (4) 

The blocking parameter a has mass dimension one. Here if), if) live on the original 
lattice while q, q live on the blocked lattice. On the second line, the integrand has been 
multiplied by one; on the last line, the if), -^-integration has been carried out. The 
induced operator on the coarse lattice is D\, and G{ x = D + aQ^Q. We will assume 
that the blocking step replaces each 2 4 hypercube by one site of the coarse lattice. 
The blocking kernel corresponds to the arithmetic mean of if) over a 2 4 hypercube, 

(QiP)(x) = 2~ A Yl V>(2av + r„). (5) 

»>=0,1 

If we apply the blocking transformation n times, the result is similar to Eq. (J3J) with 
G\ — > G n , Di — > D n , where (see ref. [TU] for further explanations) 

G- 1 = D + a n QlQ n , D- 1 = a" 1 + Q n D^Ql . (6) 

Here Q n = Q(«)g(™-i) . . . QW an d l/a n = 1/aW + l^ja^-^ + 2- 8 /a( n " 2 ) + ■■■, 
where and are respectively the blocking kernel and parameter used in the 
j-th step. 

In order to establish the locality of the fourth root of the free staggered operator 
we begin with the taste representation ([2)1 . The lattice spacing is 2a. Applying n 
blocking transformations, we obtain a coarse-lattice spacing a c = 2 n+1 a. The blocking 
parameter is fixed to be a^' = a where a = 0(1/ a c ), for all steps. We will take the 
limit n — > oo while holding a c fixed. Thus the original lattice spacing vanishes like 
2-(n+i)_ p or an y fixiite n, the generalization of Eq. (@J) gives 

det(D ) = det(G- 1 )det( J D n ). (7) 

It can be shown |6J that D n , G n , and G~ l , are all local operators on the coarse lattice, 
i.e. their kernels decay exponentially with an 0(l/a c ) decay rate, uniformly in n. The 
limit n — > oo exists, and we find det(_D ) — detf^G^ 1 ) det(D 0O ), in obvious notation. 
In the massless case the blocked propagators have the general form 

D-\p) = a' 1 - J2 fa ® QAtfp) + [ 7 5 ® nr^p)) , (8) 

where a n = 15a/(16(l — 2 -4 ™)). For small pa c , one has A™(p) = (p M /p 2 )(l + 0{pa c )). 
This corresponds to a wave-function renormalization equal to unity, as it must be 
for a free theory. In units of a c , the taste- violating amplitudes B™(p) scale like 
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a/a c oc 2~ n , uniformly in p. The taste breaking originates from an irrelevant op- 
erator (the last term in Eq. ({%])). and this is indeed the scaling anticipated for such 
an operator. Hence B^{p) = 0, and factorizes as = D Tg <g> I. The fourth root 
is det 1//4 (-Doo) = det(_D rg ). Being a fixed-point operator (in the massless case), D Tg 
satisfies the Ginsparg- Wilson relation [TT] . 

The operator G^ 1 has an 0(l/a c ) gap. The Dirac operator D has a zero at 
p = (only) and, for any n, this zero is lifted when adding the blocking-kernel part 
whose maximum is obtained at p = (see Eq. (JBJ)). As a result, there exists a local 
fourth-root operator M. satisfying det 1//4 (G^ 1 ) = det(At) having a similar gap, whose 
kernel has an 0(l/a c ) decay rate too. 

3. Interacting staggered fermions. Writing the interacting, one-component staggered 
action as x{D s + ^)x? the fourth-root trick amounts to the prescription (keeping to 
a one-flavor theory for notational simplicity) 

(X(x)x(y)---) = Z- 1 JvUdet 1/A (D s + m) e^ {D s +m)~ l {x, y) ■ ■ ■ , (9) 
Z = J VUdet 1/4 (D s + m) e~ Sg . (10) 

An essential element of the RG program is the emergence of the intermediate 
scale a c , which defines the separation between what we consider as short-distance 
and long-distance physics. In the interacting theory we will again aim at generating 
through RG blocking a new lattice theory, whose spacing a c satisfies a<Ca c < ^qcd- 
The limit n —>■ oo at fixed a c (where a/a c — > 0) will again, hopefully, give rise to a 
local lattice theory on the coarse lattice that, at the same time, will reproduce the 
physical observables of the original rooted staggered theory. We now list the main 
issues that arise in the interacting theory. A detailed account will appear elsewhere 

[TJ. 

• Gauge-field blocking and gauge invariance. The gauge-field blocking can be done in 
the usual way, see e.g. ref. The blocking of all fields must respect gauge invariance. 
For example, suitable parallel transporters must be introduced into the fermion 
blocking kernel. 

• Hypercubic invariance. For blocking of 2 4 hypercubes, it is not possible to construct 
gauge-covariant blocking kernels that will transform covariantly under hypercubic 
rotations as well. A projection onto hypercubic invariant observables can nevertheless 
be enforced. It is equivalent to introducing at each blocking step an action-less local 
field that lives on the coarse lattice; its (discrete) value at a given hypercube sets the 
relative coordinates of the site to which we will parallel transport all the variables 
residing inside that hypercube. 

• Observables. The observables of the coarse-lattice theory form a subset of the origi- 
nal observable. The RG transformation(s) define a local mapping of the coarse-lattice 
observables back into the original fine lattice. This relation is important because, in 
the absence of a normal path integral, the fourth-root prescription is defined from the 
outset in terms of observables, cf. Eq. (J§J) above. 

• From one- component to taste. As already mentioned, the interacting theory is 
originally defined in the one-component formalism. We use the resemblance between 
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Eqs. (JIJ and (0) to perform a first, special, blocking step. This step realizes a covariant 
version of Eq. (jT|). thus defining a taste-basis version of the interacting theory fUJ 
IT3*] . Since the observables of this interacting theory form a subset of the original 
observables, in effect all the symmetries of the one-component theory remain intact. 
The disastrous mass terms, which are generated if one gauges the taste representation 
in a simple-minded way [T3|, are avoided. 

• "Postponing" multi-fermion interactions. Equation (J2J) plays the key role of pro- 
viding the short- distance - long-distance separation at the level of the fermion deter- 
minant. A technical difficulty is that, as soon as we RG block the gauge fields even 
once, multi-fermion interactions are generated, and any fermionic path integral is no 
longer a determinant. The simple solution, valid for any finite number of blocking 
steps, is to perform all the fermion blocking steps ahead of all the gauge-field blocking 
steps. In this context one can still make use of Eq. (J7j). 

• Scaling. It will be unrealistic to expect rigorous proofs in the interacting case. 
With all the above elements, what one achieves is a framework for constructing an 
RG blocked version of the interacting, rooted staggered theory. The familiar scal- 
ing behavior, derived using perturbation theory, should then apply for large-enough 
number of blocking steps n, (equivalently, small-enough a). This should imply in 
particular that all the taste-breaking effects, that always originate from irrelevant 
operators, will again be damped like a/a c oc 2~ n . While not a proof, this suggests 
that, for n — > oo, the RG blocked lattice theory emerging from the rooted staggered 
theory is local. 

4. Concluding remarks. An advantage of the RG program is that, since it is formu- 
lated in a completely non-perturbative language, it is amenable to numerical tests. 
Numerical evidence for the (approximate) four-fold degeneracy of the low-lying stag- 
gered eigenmodes may be found in ref. ^Hj- For first numerical tests of the scaling 
properties of the RG blocked fermion propagator in the interacting rooted theory, see 
ref. |T3]. 
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